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O ■ ABSTRACT 

CZ3 , Context. Accretion and outflow processes are of fundamental importance for our understanding of the formation of stars and planetary systems. 

I ^1 ■ To trace these processes, diagnostic spectral lines such as the Bry 2. 166 yum line are widely used, although due to a lack of spatial resolution, 
the origin of the line emission is still unclear. 

^-H ■ Aims. Employing the AU-scale spatial resolution which can be achieved with infrared long-baseline interferometry, we aim to distinguish 
^ ' between theoretical models which associate the Bry line emission with mass infall (magnetospheric accretion, gaseous inner disks) or mass 

0^ , outflow processes (stellar winds, X-winds, or disk winds). 

^~~^ Methods. Using the VLII/ AMBER instrument, we spatially and spectrally (/i/A/l = 1500) resolved the inner (< 5 AU) environment of five 

'—' Herbig Ae/Be stars (HD 163296, HD 104237, HD 98922, MWC297, V921 Sco) in the Bry emission line as well as in the adjacent contin- 
, , uum. From the measured wavelength-dependent visibilities, we derive the characteristic size of the continuum and Bry line-emitting region. 

t ' Additional information is provided by the closure phase, which we could measure both in the continuum wavelength regime (for four objects) 

^^ ' as well as in the spectrally resolved Bry emission line (for one object). The spectro-interferometric data is supplemented by archival and new 

00 VLT/ISAAC spectroscopy 

^— ' Results. For all objects (except MWC 297), we measure an increase of visibility within the Bry emission line, indicating that the Bry-emitting 
^ , region in these objects is more compact than the dust sublimation radius. For HD 98922, our quantitative analysis reveals that the line-emitting 

• '^ ' region is compact enough to be consistent with the magnetospheric accretion scenario. For HD 163296, HD 104237, MWC 297, and V921 Sco 

Itn] ' we identify an extended stellar wind or a disk wind as the most likely line-emitting mechanism. Since the stars in our sample cover a wide range 
of stellar parameters, we also search for general trends and find that the size of the Bry-emitting region does not seem to depend on the basic 
stellar parameters (such as the stellar luminosity), but correlates with spectroscopic properties, in particular with the Ha line profile shape. 
Conclusions. By performing the first high-resolution spectro-interferometric survey on Herbig Ae/Be stars, we find evidence for at least 
two distinct Bry line-formation mechanisms. Most significant, stars with a P-Cygni Ho- line profile and a high mass-accretion rate seem to 
show particularly compact Bry-emitting regions (-^Bry/^cont < 0.2), while stars with a double-peaked or single-peaked Ha-line profile show a 
significantly more extended Bry-emitting region (0.6 < RBiy/Rmm :£ 1-4), possibly tracing a stellar wind or a disk wind. 

Key words, stars: pre-main-sequence - stars: winds, outflows - stars: individual: HD 163296, HD 104237, HD 98922, MWC 297, V921 Sco - 
planetary systems: protoplanetary disks - accretion, accretion disks - Line: formation - techniques: interferometric 

1. Introduction 



U 



Send offprint requests to: skraus@mpifr-bonn.mpg.de Accretion disks around young stellar objects (YSOs) are at the 

* Based on observations made with ESO telescopes at the La Silla focus of astronomical research, not only because they play an 

Paranal Observatory under open time programme IDs 077.C-0694, essential role in the star-formation process, but also because 

078.C-0360, and 078.C-0680. they provide the stage where planet formation takes place. 
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Historically, these disks were discovered due to their char- 
acteristic infrared excess emission, which is believed to trace 
dust grains, providing (for the expected temperatures and den- 
sities) the dominant source of opacity. In recent years, there 
has been substantial progress in constraining the detailed three- 
dimensional geometry of the dust disk using, for example, com- 
bined modeling of the spectral energy distribution (SED) and 
spatially resolved infrared interferometric observations. While 
the thermal emission from the dust disk is likely to be the dom- 
inant contributor to the infrared excess emission observed to- 
wards YSOs, it is believed that the dust content makes up only 
a small fraction of the total mass of the system. 99% of the 
mass is likely contributed by gas, in particular hydrogen, and 
can mainly be traced by the spectral lines, which for pre-main- 
sequence stars are often found in emission. While some infor- 
mation about the kinematics of the gas can already be extracted 
from the line profile, the spatial origin of the gas emission and 
the physical processes they trace are still strongly debated. 

In this context, the Brackett-y (Bry) 2.1661 yum line is of 
special importance. It was found that the luminosity L(Bry) 
of this line (as determined from the circumstellar compo- 
nent of the line equivalent width) seems to correlate with the 
mass accretion luminosity Lace, as determined from UV veil- 
ing. This empirical L(Br7)-Lacc relationship has been estab- 
lished for pre-main-sequence stars of various masses, includ- 
ing brown dwarfs (Natta et al. 2004), T Tauri stars (Muzerolle 
et al. 1998b), as well as intermediate-mass Herbig Ae/Be stars 
(Calvet et al. 2004; van den Ancker 2005). Based on this empir- 
ical correlation, observers started to use the Bry luminosity as 
an estimator for the mass accretion rate (e.g. Garcia Lopez et al. 
2006), which follows from Lace using Mace = LdccR*/GMi,. 

Therefore, it is now essential to identify the process(es) in- 
volved in the formation of the Bry emission line in YSOs. The 
main scenarios which have been proposed up to now include 
mass infall, as well as mass outflow mechanisms: 

(a) Magnetospheric accretion: The line emission might 
emerge from matter which is accreted onto the star through 
magnetospheric accretion columns (van den Ancker 2005). 
This infall is supposed to happen very close to the star, in- 
side the co-rotation radius, where the Keplerian angular veloc- 
ity matches the stellar angular velocity. 

(b) Gaseous inner disk: Inside of the dust destruction ra- 
dius, the gas continues to accrete towards the star, forming 
a gaseous inner accretion disk. The recombination line emis- 
sion from ionized hydrogen in this disk might contribute to the 
line emission observed towards HAeBe stars. Muzerolle et al. 
(2004) estimated that especially for high mass-accretion rates 
Mace > 10"^ M0yr"', the flux contribution of the inner gaseous 
disk might be of importance. 

(c) Stellar wind: The P Cygni line profile observed in the hy- 
drogen lines of several HAeBe stars might also indicate mass 
loss through stellar winds (e.g. Mestel 1968; Catala & Kunasz 
1987; Strafefla et al. 1998). While magnetically accelerated 
stellar winds seem plausible for massive stars rotating close to 
their break-up velocity, this scenario does not seem to work for 
lower-mass YSOs (Ferreira et al. 2006). 

(d) Stellar-field driven wind (X-wind): In the X-wind model 
(Shu et al. 1994), the outflows from YSOs are launched from 



a narrow interaction region where the stellar magnetosphere 
truncates the accretion disk. Keplerian rotation leads to a wind- 
ing up of the field lines and to the formation of magnetic sur- 
faces. Charged particles from the stellar wind will then get 
trapped in these surfaces, accelerated, and collimated into a 
beam. 

(e) Disk-field driven wind (disk wind): Another mechanism 
which has been proposed for the launching and collimation 
of outflows and jets observed towards YSOs are magneto- 
centrifugally driven disk winds (Blandford & Payne 1982; 
Pudritz & Norman 1983; Ferreira 1997). Material from the 
disk surface of the rotating disk is centrifugally accelerated 
along the open magnetic field lines of the accretion disk. In 
contrast to the X-wind model, where the wind is launched from 
a narrow disk annulus around the inner truncation radius of the 
accretion disk, disk winds originate from a wide range of radii, 
likely extending from the co-rotation radius out to several AUs. 

Since most of these processes are believed to take place 
on (sub-)AU scales and were therefore not accessible with di- 
rect imaging techniques, most earlier studies tried to constrain 
the spatial distribution and kinematics of the line-emitting gas 
from the shape of line profiles. However, since these line pro- 
file fitting techniques are known to be highly ambiguous (e.g. 
Catala et al. 1999), spatially resolved observations are urgently 
required to constrain the physical mechanism and spatial origin 
of important tracer lines like Bry. Since the different emitting 
mechanisms noted above can be associated to distinct spatial 
regions (see Fig. 1), spatially resolved observations should al- 
low us to identify and physically characterize the true underly- 
ing emitting mechanism. 

The VLTI/AMBER instrument combines, for the first time, 
the milli-arcsecond spatial resolution achievable with IR long- 
baseline interferometry with a good spectroscopic resolution 
(R = A/AA = 1500 or 12000) in the near-infrared /T-band. In 
this study, we use the unique capabilities of this instrument in 
order to measure the geometry and position of the Bry line- 
emitting region relative to the continuum-emitting region, pro- 
viding direct information about the processes involved. 

A particularly well-studied class of YSOs are the Herbig 
Ae/Be stars (HAeBes). These are intermediate-mass, pre-main 
sequence stars. Over the last few years, the continuum emis- 
sion from a rather large number of HAeBes could already be 
studied with broad-band infrared interferometry. These studies 
revealed a correlation between the size of the NIR continuum 
emitting region and the stellar luminosity (Monnier & Millan- 
Gabet 2002), suggesting that for most Herbig Ae stars, the NIR 
emission likely traces dust at the dust sublimation radius. 

Only very little is known about the spatial origin of the 
hydrogen recombination line emission. First AMBER obser- 
vations of the Herbig Be star MWC 297 have shown that the 
Bry-emitting region around this star is more extended than 
the continuum region (Malbet et al. 2007). Surprisingly, ob- 
servations on the less luminous Herbig Ae star HD 104237 
did not show any change in visibility along the Bry line 
(Tatulli et al. 2007a). Eisner (2007) found a visibility increase 
within the Bry line for the Herbig Ae star MWC 480. This 
might suggest that the Bry line traces fundamentally diff'er- 
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Fig. 1. Illustration of the regions which have been proposed as the origin of the permitted hydrogen recombination line emission observed 
towards HAeBe stars (this sketch is not to scale; read Sect. 1 for details about the individual mechanisms). 



ent mechanisms for Herbig Ae and Be stars. In this study, 
we present observations on three HAeBes, which are resolved 
for the first time with spectro-interferometry (HD 163296, 
HD 98922, V921 Sco). For HD 104237, we present new obser- 
vations which were obtained at longer baseline lengths, adding 
substantial new information to the initial AMBER results pre- 
sented by Tatulli et al. (2007a). Finally, we re-reduced the 
archival data on MWC 297 and HD 104237 in order to provide 
a uniformly calibrated sample of five stars for our interpreta- 
tion. The stellar parameters which we adopted for the stars in 
our sample are shown in Tab. 1 . 

This paper is structured as follows: In Sect. 2 we discuss 
the spectro-interferometric and spectroscopic data and the ap- 
plied data reduction procedures, followed by an outline of the 
obtained results (Sect. 3). Then, we present the models which 
were fitted to the data in order to constrain the geometry of the 
continuum- (Sect. 4) and Bry-emitting region (Sect. 5). After 
interpreting the individual objects in detail (Sect. 6), we discuss 
general trends which we find in our small sample of HAeBe 
stars (Sect. 7). Finally, we conclude with a summary of our re- 
sults and discuss the potential of future observations in Sect. 8. 

2. Observations and data reduction 

2.1. VLT/ISAAC spectroscopy 

The profile of the emission lines carries important information 
about the kinematics of the emitting gas. Therefore, we com- 
plement the spatially resolved AMBER spectro-interferometry 
with high-spectral resolution (R ~ 9000) spectra obtained with 
the VLT/ISAAC instrument. 

Besides archival ISAAC data (ESO programme 073. C- 
0184, RI. E. Habart), we also obtained new spectroscopic data 
(for HD 163296 and V921 Sco, ESO programme 077.C-0694, 
P.I. S. Kraus) in order to measure the line profile as close in 
time to the AMBER observations as possible. This new spec- 
troscopic data covers not only the Bry 2.1661 jjm line (Fig. 2, 
left), but also the Pa^ 1.2822 yum line (Fig. 2, right). 

The raw spectra were extracted using IRAF procedures 
and then corrected for atmospheric features using telluric stan- 



Table 3. Calibrator star information for the interferometric observa- 
tions presented in Tab. 2. 



Star 


V 


K 


Spectral 
Type 


dvD 

[mas] 


HD 101328 


7.44 


3.75 


K4III 


1.00 ± 0.01 <"' 


HD 118934 


7.92 


4.02 


K4III 


0.89 ± 0.01 <"> 


HD 135382 


2.88 


2.53 


AlV 


1.10 ± 0.08 <*' 


HD 159941 


7.85 


3.55 


MOIII 


1.09 ± 0.02 <"' 


HD 171960 


7.29 


3.38 


K3II 


1.13 ± 0.02 <"> 


HD 177756 


3.43 


3.56 


B9Vn 


0.60 ± 0.06 **> 



Notes - The V-band magnitudes were taken from SIMBAD and 
the A'-band magnitudes from the 2MASS point source catalog. 

References - (a) UD diameter taken from the CHARM2 
catalog (Richichi et al. 2005). (b) UD diameter computed with 
ASPRO. 



dard star observations obtained during the same night. For 
wavelength-calibration, we aligned the raw spectra to publicly 
available high-resolution (R - 40000) telluric spectra taken at 
the NSO/Kitt Peak Observatory. To compare the final ISAAC 
spectra (Fig. 2) with the spectra extracted from the AMBER 
data, we convolved the ISAAC spectra to the spectral resolu- 
tion of AMBER (R - 1500) and found good agreement (see 
Figs. 3(a) to 7(a), top panel). 



2.2. Archival ISO and Spitzer/IRS spectroscopy 

In order to optimally constrain the SED for the stars in our sam- 
ple, we obtained near- to mid-infrared spectroscopic data from 
the ISO and Spitzer Space Telescope archive. The Spitzer/IRS 
spectra were pre-processed by the S13.2.0 pipeline version at 
the Spitzer Science Center (SSC) and then extracted with the 
SMART software, Version 6.2.5 (Higdon et al. 2004). 
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Table 1. Target stars and adopted stellar parameters. 



Star 



Spectral 
Type 






d 
[pc] 



T, 
[K] 



Av 



[Rq} 



log ii222 (ft) 



logM^J*^' 
[Mgyr-'l 



V sin I 
[kms-'] 



Ha'* 
profile 



Ref. 



HD 104237 


A5 


30 


116 <^> 


8 000 


0.31 


2.9 


-2.74 


-7.45 


12 ± 2 <s> 


DC-) 


(0 


HD 163296 


A3 


26 


122 (/) 


8 700 


0.12 


2 2 


-2.78 


-7.12 


120+20 U) 


D(*) 


(k) 


HD 98922 


B9 


890 


540 <'^' 


10600 


0.3 


9.1 


-1.56 


-5.76 


- 


p(A) 


(0 


MWC 297 


B1.5 


10600 


250 


23 700 


8 


6.1 


>-0.6 


- 


350 + 50 <""' 


S<*' 


(m) 


V921 Sco 


BO 


19950 


800 


30900 


5.0 


5.0 


>-0.7 


- 


- 


D<"' 


(«) 



Notes - *"' The stellar radius Rf, is computed using the given eifective temperature and bolometric luminosity. 

** Luminosity of the Bry emission line, as determined by Garcia Lopez et al. (2006) from ISAAC spectra. 

*''' Mass accretion rate, as determined from the L(Br7)-Lacc relation (Garcia Lopez et al. 2006). 

** This column described the Ha-line profile shape of the stars in our sample using the classifications: D: double-peaked profile; P: 
P-Cygni profile; S: single-peaked profile. 

*'"' For HD 98922, we assume the minimum distance given by Garcia Lopez et al. (2006). 

References: (/) van den Ancker et al. 1998; (g) Donati et al. 1997; {h) Acke et al. 2005; (0 TatuUi et al. 2007a; (j) Finkenzeller 1985; 
ik) van den Ancker et al. 2000; (/) Garcia Lopez et al. 2006; (m) Drew et al. 1997; {n) Habart et al. 2003 



Table 2. Observation log of our VLTI/ AMBER interferometric observations with a spectral resolution R - 1500. 



Target Star 



Date 

(UT) 



Time 

(UT) 



Spectral 

Window 

[}im\ 



DIT 

[ms] 



Telescope 
Triplet 



Projected Baselines 



[m] 



PA 



Bl 
[m] 



PA 

n 



Bl 

[m] 



PA 

n 



Calibrator Ref. 



HD 104237 


2005-02-26 


07:14 


2.12-2.20 


100 


HD 104237 


2007-01-09 


08:03 


1.94-2.26 


200 


HD 104237 


2007-01-09 


07:41 


1.94-2.26 


500 


HD 163296 


2006-04-13 


06:19 


2.12-2.20 


50 


HD 98922 


2007-02-04 


07:19 


2.12-2.20 


50 


HD 98922 


2007-01-09 


09:31 


1.94-2.26 


200 


HD 98922 


2007-01-09 


09:39 


1.94-2.26 


500 


MWC 297 


2004-05-31 


06:02 


2.00-2.23 


107 


V921 Sco 


2006-04-14 


08:09 


2.12-2.20 


50 



UT2-UT3-UT4 
UT1-UT3-UT4 
UT1-UT3-UT4 
UT2-UT3-UT4 

UT2-UT3-UT4 
UT1-UT3-UT4 
UT1-UT3-UT4 

UT2-UT3 
UT2-UT3-UT4 



35.1 


71 


(78.0 


31) 


(78.7 


26) 


43.9 


19 


42.9 


46 


(90.5 


42) 


(89.9 


43) 


44.7 


42 


45.2 


42 



(60.8 


120) 


58.5 


84 


58.0 


78 


52.3 


99 


(61.9 


111) 


62.3 


116 


62.4 


118 


62.1 


110 



(87.9 


102) 


HD 135382 


( 122.4 


53) 


HD 118934 


(123.1 


47) 


HD 118934 


74.1 


63 


HD 171960 
HD 89682 


(88.9 


85) 


HD 101328 


(122.8 


71) 


HD 118934 


(122.0 


73) 


HD 118934 


- 


- 


HD 177756 


(89.3 


82) 


HD 159941 



(a) 



(b) 



Notes - For more detailed information about the calibrator stars, we refer to Tab. 3. 

References - (a) reprocessing of data presented in Tatulli et al. (2007a); (b) reprocessing of data presented in Malbet et al. (2007). 



2.3. VLTI/AMBER spectro-interferometry 

AMBER (Petrov et al. 2007) is the NIR beam-combiner of 
the Very Large Telescope Interferometer (VLTI), which is lo- 
cated on Cerro Paranal/Chile and operated by the European 
Southern Observatory (ESO). Combining the light from up 
to three of the four 8.4 m unit telescopes simultaneously, 
AMBER measures not only visibility amplitudes, but also 
the closure phase (CP) relation. In the course of three ESO 
open time programmes (077.C-0694, 078.C-0360, 078.C-0680, 
P.I. S. Kraus), we obtained spectrally dispersed interferograms 
in AMBER'S medium resolution (MR) mode (R - 1500) 
on four HAeBes; namely, HD 163296, HD 104237, HD 98922, 
and V921 Sco. These new data sets were complemented with 
archival data from MWC 297 and HD 104237, obtained earlier 
during AMBER commissioning and in guaranteed-time obser- 
vations (Malbet et al. 2007; Tatulli et al. 2007a). Re-reducing 
these data sets with the latest software allows us to ensure ho- 
mogenity both in data reduction as well as in the applied mod- 



eling procedures. The observations are summarized in Tab. 2 
and were obtained under good atmospheric conditions (seeing 
0.5-0.9", atmospheric coherence time 2-7 ms). In order to cal- 
ibrate the obtained visibilities and CPs for instrumental and 
atmospheric effects, we used interferometric calibrator stars, 
taking the intrinsic diameters of these stars into account (see 
Tab. 3). By using different detector integration times (DITs) 
for the different objects, we aimed for a compromise between 
collecting a sufficient number of photons and minimizing the 
loss of fringe contrast due to atmospheric piston. In general, 
observations with longer DIT should provide better SNR in 
the differential visibility measurement, while observations with 
shorter DIT should provide a better absolute calibration. Since 
we aimed mainly for a precise measurement of the differen- 
tial quantities in the Bry and the adjacent continuum, we used 
rather long DIT for most observations. In order to account for 
the resulting errors on the absolute calibration, we estimate a 
rather large calibration error of 5%. The spectral calibration of 
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the data was done by comparing the spectrum extracted from 
the AMBER data with ISAAC spectra of the same object (see 
Sect. 2.1), which have been convolved to the spectral resolution 
of AMBER. An additional very important calibration step in 
the course of AMBER data reduction was the relative spectral 
re-shifting of the photometric channels with respect to the in- 
terferometric channel. The appropriate shift was determined to 
sub-pixel accuracy by computing the auto-correlation between 
the spectra extracted from these channels. 

In a first data reduction step, the AMBER raw interfero- 
grams were cleared from the correlated detector noise effect 
using the AMDC software tool (Version 1.1, Causi et al. 2007). 
Then, the raw data was reduced with the amdlib2 software^ 
(release 2.1), employing the P2VM algorithm (Tatulli et al. 
2007b). Due to the absence of a fringe tracker, a large fraction 
of the interferograms is of rather low contrast (see discussion 
in Petrov et al. 2007). Therefore, we removed any frames from 
our dataset for which the light injection from the contributing 
telescopes was unsatisfying; i.e., the intensity ratio between the 
photometric channels was larger than 4, or the fringe contrast 
was decreased due to instrumental jitter (the 20% best inter- 
ferograms were selected based on the Fringe SNR criteria, as 
defined in TatulU et al. 2007b). 

Since some stars in our sample are close to the current sen- 
sitivity limit of amber's MR mode, we applied spectral bin- 
ning to the raw data in order to increase the fringe SNR. For 
each data set, the width of the sliding window was chosen so 
that the resulting fringe SNR of the interferograms exceeds the 
critical value of 1.5, ensuring a reliable visibility estimation. 
Since the Bry line for all our objects is spectrally resolved over 
several spectral channels ai R - 1500, the resulting decrease 
in spectral resolution only marginally reduces the contrast be- 
tween the spectral line and the underlying continuum. This pro- 
cedure increases the fringe SNR significantly, while the loss 
in spectral resolution (from R - 1500 to 500 or 250) results 
only in a minor decrease of the line-to-continuum flux ratio 
^Biy/^cont- BascHnes for which a spectral binning io R - 250 
was not sufficient to yield an SNR of 1.5 were rejected from 
further analysis (in Tab. 2, these baselines are put in brackets). 
The final visibility curves are shown in Figs. 3(a) to 1(a) in the 
second panel from the top. 

In principle, two distinct phase quantities can be extracted 
from AMBER interferograms, namely differential phases (DP, 
measuring the relative displacement of the line-emitting region 
with respect to the continuum-emitting region) and the closure 
phase, which can indicate asymmetries in the source bright- 
ness distribution. Since for most of our datasets, the fringe 
SNR in individual frames is not sufficient to reliably correct 
for disturbing atmospheric phase contributions (atmospheric 
piston), no scientifically meaningful differential phases could 
be extracted, unfortunately. However, for some data sets we 
could extract useful CP information, both in the continuum 
(HD 104237, HD 163296, HD 98922, V921 Sco) as well as in 
the line regime (V921 Sco). It is expected that the continuum 
CP does not change significantly over the small spectral win- 



dow covered by our observations, which allowed us to apply an 
optimized data reduction strategy for the extraction of the con- 
tinuum CPs. First, we apply some spectral binning to the raw 
data (bracketing out the spectral channels containing the emis- 
sion line), followed by an averaging of the complex visibilities 
over various spectral channels (between 2.12-2.15 and 2.175- 
2.185 jum). The obtained CPs will be presented in Sect. 3.3. 

3. Results 

3.1. SED and NIR-spectroscopy 



velocity [km/s] 

■250 250 500 



velocity [km/s] 

■500 -250 250 




' The amdlib2 software package is available from the website 
http : //www . jmmc . fr/data_processing_amber . htm 



2.162 2.164 2.166 2.168 2.17 

Wavelength [|xm] 

Fig. 2. Left: VLT/ISAAC spectra showing the Bry line with a spectral 
resolution ofR ~ 9000 for the stars in our sample. Right: VLT/ISAAC 
spectra of the Pa/3 line were obtained for HD 163296 and V921 Sco. 



In order to constrain the SED of the stars in our sample, we 
collected photometric and spectroscopic data from the litera- 
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ture and the ISO and Spitzer archives (see Sect. 2.2). Assuming 
that the line-of-sight extinction can be mainly accounted to the 
large scale envelopes, in which some HAeBe stars are embed- 
ded, we dereddened the photometric and spectroscopic data as- 
suming the values for Ay given in Tab. 1, Ry - 3.1, and the 
extinction law by Mathis (1990). The derived SEDs (Figs, iib) 
to l(bj) allow us to determine the ratio between the star flux 
contribution /star and the flux contribution from the circumstel- 
lar disk /disk for the fitting of geometric models (Sect. 4). 

The VLT/ISAAC Bry line spectra are shown in Fig. 2. With 
the currently available resolution, all Bry lines show a single- 
peaked profile. For HD 163296 and V921 Sco, measurements 
at two or three epochs are available, indicating that the line 
profile of V921 Sco has not significantly changed, while for 
HD 163296 we find some variability. 

3.2. Spectro-interferometric visibilities 

The wavelength-dependent visibilities extracted from our 
AMBER data are shown in the second panel from the top 
in Figs. 3(a) to 7(a). Four stars in our sample (HD 163296, 
HD 104237, HD 98922, V921 Sco) show a visibility increase 
within the Bry line. For MWC 297, the visibility drops within 
the line, as already reported by Malbet et al. (2007). 

Since the spectral channels which include the Bry line 
additionally contain flux contributions from the photosphere 
and circumstellar material, the interpretation of the measured 
wavelength-dependent visibilities requires quantitative model- 
ing, as presented in Sect. 4. 

3.3. Closure phases 

Closure phase measurements can provide unique information 
about deviations from centro-symmetry in the source bright- 
ness distribution. For YSO disk geometries, such asymmetries 
are expected in particular for systems seen under intermedi- 
ate inclination. The strongest CP signals are predicted by disk 
models with vertical puffed-up inner rim (e.g. Dullemond 
et al. 2001), while models with curved inner rims predict a 
smoother, more symmetric brightness distribution correspond- 
ing to smafler CPs (e.g. Isefla & Natta 2005). 

Given the importance of this observable, we present here 
the CPs extracted from our data, although the error bars on 
most measurements are rather large, typically due to the low 
fringe contrast on the longest baseline in the telescope triplet. 
In Tab. 4 we list the CPs measured from our datasets and 
give the corresponding (m, v) baseline vector for two base- 
lines in the employed telescope triplet (the third (m, v)-vector 
is given by the closure relation). We find that all continuum 
CPs are consistent with a zero CP on the Icr (HD 163296, 
HD 98922, V921 Sco) or 2cr level (HD 104237). This adds sup- 
port to the conclusions drawn by Monnier et al. (2006), who 
used the I0TA-3T interferometer (providing baseline lengths 
up to 38 m) and measured for 12 out of 13 sources CPs be- 
low ~ 5° (excluding one binary source). Since the CP signal is 
expected to increase rapidly with baseline lengths, our new ob- 
servations (with baselines up to 89 m) provide new constraints 



on this issue and favour again a centro-symmetric brightness 
distribution. In Sect. 6, we discuss some of our CP measure- 
ments qualitatively, but leave it to future studies to investigate 
whether our results are also in quantitative agreement with the 
current generation of rim models. 

For V921 Sco, the particularly strong Bry line flux al- 
lowed us also to measure an accurate CP measurement within 
the spectrally resolved emission line, yielding a CP signal of 
3.1 ±2.8°, consistent again with a centro-symmetric brightness 
distribution. 

4. Modeling: Continuum-emission 

Over the last few years, ring models with uniform brightness 
have emerged as the prototypical geometry for the interpreta- 
tion of YSO interferometric data. In most cases, this prefered 
use of ring geometries instead of other simple geometries (e.g. 
uniform disk, Gaussian, or disk geometries with constant tem- 
perature power-law) is not explicitly required to reproduce the 
details of the sampled visibility function (which typically cov- 
ers only the first lobe and is thus rather insensitive to the inner 
gap in the brightness distribution), but mainly based on indi- 
rect evidence or theoretical arguments. For example, some of 
the first infrared interferometric measurements on HAeBe stars 
(e.g. Millan-Gabet et al. 2001) already showed that classical 
geometrically thin accretion disk geometries extending down 
to several stellar radii might be consistent with interferometric 
measurements for early-type Herbig Be stars, but result in too 
compact structures for Herbig Ae and late-type Herbig Be stars 
(Eisner et al. 2004; Vinkovic & Jurkic 2007). Ring geometries, 
on the other hand, yield sizes which are consistent with the 
expected dust sublimation radii and follow the predicted stel- 
lar luminosity scaling law (Monnier & Millan-Gabet 2002). In 
the meantime, this finding was also interpreted in theoretical 
work (e.g. Natta et al. 2001; Dullemond et al. 2001; IseUa & 
Natta 2005), attributing the near-infrared continuum emission 
mainly to hot dust located at the dust sublimation radius. Based 
on these arguments, we also prefer to use ring-like geometries 
for the interpretation of our visibility data (assuming a frac- 
tional ring width of 20%, Monnier et al. 2005), but also give 
uniform disk (UD) and Gaussian FWHM diameters to allow 
comparison with other work. In order to estimate the contribu- 
tion of the stellar photosphere to the total TT-band flux, we es- 
timate the flux ratio using the SED and the Kurucz atmosphere 
models shown in Figs. 3(b) to 7(b). For more details about the 
model fitting procedure, we refer to one of our earlier studies 
(Kraus et al. 2008). 

Besides model fits to each individual spectral channel (see 
fourth and fifth panel of Figs. 3(a) to 7(a)), we also fitted these 
models to continuum-visibilities, which were averaged over 
several spectral channels around 2. 15 and 2.18 /im and plot the 
corresponding measurements and model curves as a function of 
baseline length (Figs. 3(c) to 7(c)). In these plots, we also show 
the upper limits, which we can put on the continuum visibility 
using the baselines which had to be rejected due to low SNR 
(see Sect. 2.3). The fitted diameters for Gaussian, UD, and ring 
profiles are listed in Tab. 5. To convert the measured angular 
size to physical scales, we assume the distances listed in Tab. 1 . 
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Table 4. Closure Phases extracted from the data. 





Projected Baselines 


Closure Phase 




"12 


V\2 


"23 


V23 


continuum 


Bry line 




[m] 


[m] 


[m] 


[m] 


[deg] 


[deg] 


HD 104237 


33.1 


11.6 


52.9 


-30.0 


9.8 + 6.7 


— 


HD 163296 


14.5 


41.4 


51.6 


-8.0 


1.5+4.3 


— 


HD 98922 


30.9 


29.7 


57.6 


-22.6 


-9± 16 


— 


V921 Sco 


30.2 


33.7 


58.3 


-21.6 


20 + 24 


3.1 ±2.8 



Table 5. Best-fit parameters for our geometric model fits. 









Continuum emission 






Bry emission 


Bry emission 
















line center 


all spectral channel 




Gaussian 




Uniform Disk 


Ring 




Ring 




Ring 






FWHM,„„ 


^red 


^cont 


^rcd 


-^cont 


2 
^rcd 


£>Bry 


2 


£»Br, 


2 


Target Star 


[AU] 




[AU] 




[AU] 




[AU] 




[AU] 




HD 104237 


0.57 ±0.10 


1.0 


0.87 ± 0.08 


0.5 


0.58 ± 0.04 


0.3 


0.31 ±0.13 


0.6 


0.35 ±0.21 


0.5 


HD 163296 


0.39 ± 0.05 


0.9 


0.61 ±0.08 


1.5 


0.41 ± 0.06 


1.9 


0.35 ±0.15 


0.1 


0.25 ±0.19 


0.1 


HD 98922 


2.45 ± 0.26 


0.5 


3.70 ± 0.27 


0.2 


2.47 ±0.16 


0.2 


<0.5 


— 


<0.5 


— 


MWC297 


1.27 ±0.10 


— 


1.99 ±0.13 


— 


1.35 ± 0.08 


— 


1.89 ±0.52 


0.1 


2.52 ± 0.44 


0.4 


V921 Sco 


5.38 ±0.52 


0.7 


8.22 ±0.55 


0.1 


5.49 ± 0.32 


0.1 


4.45 ± 0.55 


0.1 


4.35 ± 0.64 


0.2 



Notes ~x^t^i is defined as 2 [(K 



IS - V'n,odei)/o"v,„aJ /(^ ~ 1)' where A' is the number of measurements, 
are the measured visibility and total error, respectively. 



is the model visibility, and 



Please note that the errors given do not take distance uncertain- 
ties into account. By fitting the visibihties measured on dif- 
ferent basehne orientations simultaneously, we assume that the 
intensity profile does not depend on position angle (e.g. face- 
on disk). Since the visibilities derived for HD 163296 provide 
sufficient position angle coverage to investigate for a possible 
elongation of the continuum-emitting region (e.g., due to disk 
inclination), we also fitted inclined ring geometries to this data 
set and derived an inchnation of 68+10° with a PA of 144 + 9° 
(corresponding to an inclined ring diameter of 1 .02 x 0.37 AU). 

5. Modeling: Bry line emission 

From each spectral channel of our AMBER interferograms, we 
can derive a value for the visibility amplitude, providing spa- 
tial information about the brightness distribution contributing 
to this spectral channel. Within the Bry line, the flux within a 
spectral channel is composed of the line emission plus the un- 
derlying continuum contribution. To model such visibility data, 
either the composite (line & continuum) object can be modeled, 
or the measured continuum-nline visibility has to be corrected 
for the continuum contribution, yielding the visibility of the 
pure line-emitting region. 

We correct the measured visibilities using the relation 



^Bry - 



■FfntVtr 



V„ 



cont ' cont 



Bry 



(1) 



where Ftot and Vtot are the flux and the visibility measured 
within an arbitrary spectral channel. To compute the Bry- 
line flux F-^^ :- Ftot - ^cont, we determine the continuum 



flux /^cont from Kurucz atmosphere models, taking the under- 
lying photospheric Bry absorption component into account. 
Then we interpolate the continuum visibility Vcont over the 
Bry line and derive the continuum-corrected visibility Vbyj of 
the Bry-emitting region. As discussed in Weigelt et al. (2007, 
Appendix C), this relation assumes a negligible (zero) differen- 
tial phase (corresponding to coinciding photocenters between 
the continuum- and line-emitting regions). We compute the 
continuum-correction for each spectral channel of the AMBER 
data where F^^ is sufficiently large to apply Eq. 1 reliably and 
show the resulting Vfiry values in the third panel from the top 
of Figs. 3(a) to 7(a). 

The interpretation of the derived line visibilities is diffi- 
cult for several reasons. For instance, it is likely that the ge- 
ometry of the line-emitting region is more complex than the 
continuum-emitting region, possibly extending above the equa- 
torial disk plane, perhaps introducing strong inclination effects. 
Furthermore, as indicated by the short-period line variability 
detected towards YSOs of all masses, the kinematics and pos- 
sibly also morphology of the line-emitting gas might change 
significantly even on short time scales of days or weeks, in- 
volving, for example, the ejection of fast-moving blobs at the 
base of an outflow. Finally, it is likely that several processes 
(including accretion and outflow processes) contribute to the 
Bry emission. Given this complexity and the large error bars, 
we refrain from applying complex line radiative transfer mod- 
els for the interpretation of the measured line visibilities in this 
first paper and instead limit ourselves to fitting simple ring- 
geometries, assuming one dominant source for the Bry emis- 
sion. Since the derived line visibility values are rather high for 
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most objects, the characteristic size of the Bry-emitting region 
does not critically depend on the intensity profile of the un- 
derlying brightness distribution and therefore affects our con- 
clusions only marginally. Since the Bry line was spectrally re- 
solved by AMBER in all observations, we performed two fits: 
one including all spectral channels where Vsry could be deter- 
mined and one including only the spectral channel at line cen- 
ter. The circular-symmetric ring-model fits to the continuum- 
corrected Bry visibilities are shown in Figs. 3(d) to 7(d), and 
the determined best-fit ring diameters are listed in Tab. 5. 

6. Discussion of individual objects 

6.1. HD 104237: A circumbinary disk 

HD 104237 (alias DX Cha, CD-77°528, Hen 3-741) has an ap- 
proximate age of ~ 2 Myrs (Bohm et al. 2004). Imaging ob- 
servations by Grady et al. (2004) suggest that HD 104237 is 
associated with a microjet (HH 669) and that the disk is seen 
nearly face-on (/ = 18^[| °). Spectroscopic monitoring suggests 
that HD 104237 consists of a binary system on a short-period 
(19.859 d), highly eccentric (e x 0.66) orbit with component 
masses of M^ = 2.2 + 0.1 M© and Ms = 1.7 + 0.1 M© (Bohm 
et al. 2004). Using these orbital elements, we estimate the ma- 
jor axis of the binary system to be ~ 1 .91 mas (0.23 AU), which 
places both stellar components inside of the ring-diameter mea- 
sured by our AMBER observations at continuum wavelengths 
(~ 5.0 mas - 0.58 AU). The measured continuum ring radius of 
0.29 AU is very close to the expected dust sublimation radius 
(0.32 AU, L = 30 Lq, T,m = 1500 K, see Sect. 7.1), suggest- 
ing that the continuum emission traces mainly dust located at 
the inner rim of the circumbinary dust disk. 

In the Bry line, we measured no significant deviation of 
the line visibility with respect to the continuum visibility on 
2005-02-26 at a short baseline length (Bp ~ 35 m), while on 



2007-01-09, the new measurements at B„ 



58 m show a 



significant increase in visibility (see Fig. 3). In both epochs, 
the continuum-corrected Bry-line visibility indicates that the 
emitting region has a similar or smaller extension than the 
continuum-emitting region. In all available measurements, the 
Bry-emitting region is significantly more extended than the co- 
rotation radius, which leads us to reject the magnetospheric ac- 
cretion scenario and favor the disk wind scenario (consistent 
with the conclusion from Tatulli et al. 2007a). Due to the low 
mass accretion rate of this object (log Mace = -7.45, Garcia 
Lopez et al. 2006), we deem major line contributions from a 
gaseous inner disk unlikely (Muzerolle et al. 2004). 

For each epoch, the data cover only single baseline lengths, 
which makes it difficult to estimate the influence of the bi- 
nary companion both on the continuum and the Bry visibilities. 
Given our marginal detection of a non-zero continuum closure 
phase signal for HD 104237 (see Tab. 4) and the estimated bi- 
nary parameters, it seems likely that we also detect small con- 
tributions from the inner resolved binary system. In particular, 
the binary nature might result in a significant asymmetry in 
the brightness distribution (in continuum and/or line emission) 
which could bias both the correction for the stellar contribution 
(see Sect. 3.1) as well as the continuum-correction required to 



compute Vfiry (see Sect. 5 and Eq. 1). Therefore, future investi- 
gations on this object will require a significantly larger amount 
of observational data to reliably discern the influence of these 
effects. 



6.2. HD 163296 

HD 163296 (aUas MWC275) is a particularly wefl-studied 
Herbig Ae star. CO millimeter line observations have revealed 
the outer parts of a Keplerian rotating disk (Mannings & 
Sargent 1997; Isella et al. 2007). The determined disk incli- 
nation angle (i - 46 + 4°) as well as the position angle 
(PA=128 + 4 °) are in reasonable agreement with the orienta- 
tion of a disk-like structure seen in scattered light from corono- 
graphic imaging (/ = 60 ± 5°, PA=140 + 5°, Grady et al. 2000). 
Imaging observations have revealed an asymmetric, bipolar jet 
emerging roughly perpendicular to the disk plane and termi- 
nating in the Herbig-Haro objects HH 409-A to C, which were 
found along PA=42.5 + 3.5°. The inclination angle inferred 
from the jet observations is / = 51^^' ° (Wassell et al. 2006). 

The infrared continuum environment around HD 163296 
was already resolved with the VLTI/MIDI (A^-band; y(99m) ss 
0.2; Leinert et al. 2004), IOTA (//-band; y(38m) » 0.85; 
Millan-Gabet et al. 2001), and the Keck interferometer (K- 
band; y(76m) x 0.47; Monnier et al. 2005). Recent measure- 
ments obtained with the CHARA interferometer (Tannirkulam 
et al. 2008) indicate that a significant fraction of the NIR con- 
tinuum emission emerges from inside the dust sublimation ra- 
dius, maybe due to gas contributions. 

High-resolution spectroscopy (R - 25 000) of the Bry line 
of HD 163296 was presented by Brittain et al. (2007), showing 
a single-peaked line profile consistent with our ISAAC spectra. 
For the origin of the Balmer lines, Pogodin (1994) proposed a 
stellar wind/outer shell model. 

In the y(B)-plot for the continuum wavelength regime 
(Fig. 4(c)), our AMBER measurements towards different base- 
line position angles (covering a PA range of 80°) cannot be well 
represented with the visibility profiles corresponding to ring- 
, UD-, or Gaussian-geometries, which might indicate either a 
very unusual radial intensity profile or, more likely, an elon- 
gation of the continuum-emitting region. The latter hypothesis 
can be investigated with the inclined ring model fits shown in 
Fig. 4(e), indicating in fact that the near-infrared emitting re- 
gion is strongly elongated. The major axis of the continuum- 
emitting region (PA 144 + 9°) seems to be oriented perpendic- 
ular to the outflow direction indicated by the optical microjet 
(PA 42.5 + 3.5°, WasseU et al. 2006). Due to the limited po- 
sition angle coverage of our interferometric observations, the 
disk inclination angle (/ = 68°) is only poorly constrained by 
our observations, but seems to be in general in good agree- 
ment with the high inclination inferred from earlier observa- 
tions. In spite of the strong inclination, the derived continuum 
closure phase for HD 163296 is still surprisingly close to zero 
(see Tab. 4), favouring disk models without notable inclination- 
induced skews. 

As indicated by the quite flat wavelength-dependent visi- 
bility curve (Fig. 4(a), second panel), the Bry emission orig- 
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inates in a region of similar spatial extent as the continuum- 
emitting region. On the longer baselines, we find weak evi- 
dence (~ 2cr) for an increase of visibility within the Bry line, 
which indicates that most of the line emission emerges from 
within the dust-sublimation radius. Plotting the continuum- 
corrected Bry-line visibility as a function of baseline length 
shows that the Bry-emitting region is likely spatially resolved 
(with RBiy/Rcont ~ 0.7), which leads us to favor a stellar wind 
or a disk wind as the emitting mechanism. Considering the low 
mass accretion rate of HD 163296 (log Mace - -7.12, Garcia 
Lopez et al. 2006), it seems unlikely that the Bry line emission 
emerges from a gaseous inner disk (Muzerolle et al. 2004). 

6.3. HD 98922 

For the B9Ve-type star HD 98922, a moderate infrared excess 
(Malfait et al. 1998) and a very high mass accretion rate was de- 
termined (Garcia Lopez et al. 2006; log Mace - -5.76 Mq yr"', 
which is the highest mass accretion rate in the sample of 36 
Herbig Ae stars investigated by these authors). Outflow sig- 
natures were detected in the Ha as well as in the Sin line 
(v ~ 300 kms-i, Grady et al. 1996). Acke et al. (2005) tried 
to model the profile of the [O i] emission line with a Keplerian 
rotating disk model and concluded that the emission likely does 
not originate in the disk surface as in several other sources, but 
possibly in a rotating gaseous disk inside the dust-sublimation 
radius. Van Boekel et al. (2003) studied the profile of the 
10 /im silicate emission feature detected towards HD 98922 and 
found spectral features associated with olivines and pyroxenes, 
indicating the presence of large, rather evolved grains. 

The inclination under which HD 98922 is seen is still not 
well known. Blondel & Djie (2006) modeled the UV spectrum 
and derived an intermediate inclination angle (/ x 45). Since 
the ring-diameters which we derive from our AMBER mea- 
surements are nearly identical towards very different position 
angles (Fig. 5(0;), lowest panel), we find no indications for a 
large disk inclination, but cannot exclude it either. Despite the 
fact that HD 98922 has the lowest observed Bry line equiva- 
lent width of the stars in our sample {F-Bry/Fcom ~ 1-2), we 
measure a strong increase of visibility within the Bry line with 
respect to the continuum visibility (see Fig. 5). Remarkably, 
applying the continuum-correction to the measured Bry-line 
visibility reveals that the Bry line-emitting region seems to be 
completely unresolved (i.e. Vsry ~ 1) at baseline lengths up to 
~ 60 m, while the continuum-emitting region is resolved with a 
ring-diameter of 4.6 + 0. 1 mas (2.48 + 0.05 AU). Therefore, the 
Bry-emitting region is at least a factor of 5 times more com- 
pact than the dust-sublimation radius, which is consistent with 
the stellar wind, X-wind, and magnetospheric accretion sce- 
nario. Given the spectroscopic indications for an exceptionally 
high mass accretion rate, we favor the magnetospheric accre- 
tion scenario. 



6.4. MWC297 

MWC297 (alias NZ Ser) is a B1.5Ve-type Herbig Be star with 
a remarkably high rotation velocity of 350 + 50 kms"' (as 



derived from photospheric lines), suggesting that the system 
might be seen nearly edge-on (Drew et al. 1997). On the other 
hand, using spectro-polarimetry in the Ha line, Oudmaijer & 
Drew (1999) could not find any evidence for asymmetry in the 
Bry line-emitting region, as would be expected for a nearly 
edge-on system. Malbet et al. (2007) applied an analytic ac- 
cretion disk model and could reproduce the SED and AMBER 
/T-band and IOTA //-band visibilities with an inclination angle 
of / ~ 20° (i.e. closer to face-on than edge-on). On the other 
hand, Acke et al. (2008) investigated the near- and mid-infrared 
geometry of the circumstellar environment around MWC 297 
and found best agreement using a two-component Gaussian ge- 
ometric model. 

The hydrogen line emission observed towards MWC 297 is 
exceptionally strong and has been interpreted as stellar wind 
(Nisini et al. 1995; Malbet et al. 2007) or accretion-driven 
mass-loss from the disk (Corcoran & Ray 1998). 

The UKIRT spectra of Drew et al. (1997) show a single- 
peaked profile for the Ha, Hj3, and Bra line. Garcia Lopez 
et al. (2006) presented an ISAAC spectrum of MWC 297, 
showing a double-peaked Bry-line profile. Based on our re- 
reduction of the ISAAC spectroscopic data, we cannot confirm 
the double-peaked profile of the Bry line, but classify the line 
profile as single-peaked with the currently available resolution 
(R ~ 9000). NIR spectroscopy with higher spectral resolution 
will be required to determine the real underlying line profile. 

As already discussed by various other authors (Monnier 
et al. 2005; Malbet et al. 2007), the continuum visibilities 
measured on MWC 297 indicate a very compact continuum- 
emitting region, corresponding to a ring-radius approximately 
5-times more compact than the dust sublimation radius ex- 
pected for an irradiated dust disk. The interpretation of this ef- 
fect (which was also observed towards several other Herbig Be 
stars) is still debated, but might include either gas absorption 
(allowing dust to exist closer to the star; e.g. Monnier & 
Millan-Gabet 2002; Malbet et al. 2007) or the emission of an 
optically thick gaseous inner accretion disk (e.g. Monnier et al. 
2005; Kraus etal. 2008). 

For MWC 297, the visibility drops significantly within the 
Bry line (see Fig. 6), indicating that the line-emitting region 
is more extended than the continuum-emitting region. Malbet 
et al. 2007 modeled this observational result assuming an op- 
tically thick gas disk with an inner radius of 0.5 AU and an 
outflowing stellar wind. Besides this interpretation, a disk wind 
scenario would likely also be consistent with our data. In this 
context, it is important to note that the pronounced drop of vis- 
ibility within the Bry line is also supported by the particular 
compactness of the continuum-emitting region of MWC 297. 
As discussed above, the continuum-emission around MWC 297 
might not be dominated by pure dust emission, but instead by 
the emission of gas+dust located close to the star. Therefore, 
the observation that MWC 297 shows a drop of visibility within 
the Bry line, while V921 Sco (which has a similar spectral type 
as MWC 297) shows an increase of visibility, might not indi- 
cate substantially different line formation mechanisms, but may 
simply reflect different physical conditions in the continuum- 
emitting dust and gas disk (e.g. optically thick vs. optically thin 
gas disk). 
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In case of a nearly edge-on system inclination, the influ- 
ence of geometric effects (extension of the Bry-emission region 
above the mid-plane, appearance of bipolar lobes, etc.) would 
also be of importance for the interpretation and modeling of the 
spectro-interferometric data on MWC 297. 



6.5. V921 SCO 

In the catalog of HAeBe member and candidate stars by 
de Winter & The (1990), V921 Sco (ahas MWC 865, CD- 
42°11721, Hen 3-1300) is classified as a BO[e]p-type star, al- 
though there is some dispute both about the stellar parame- 
ters and the evolutionary stage of this object (see discussion 
in Borges Fernandes et al. 2007). In the vicinity of V921 Sco, a 
cluster of YSO candidate sources were detected (Habart et al. 
2003; Wang & Looney 2007), and the star is also associated 
with an extended nebulosity, which was observed at optical 
(van den Bergh & Herbst 1975; Hutsemekers & van Drom 
1990), near-infrared (Wang & Looney 2007), mid-infrared 
(Natta et al. 1993), as well as sub-mm wavelengths (Henning 
et al. 1998). Comparing the spectral slope of the ISO spectrum 
with that of the Spitzer/IRS spectrum taken with a smaller aper- 
ture (see Fig. 7b) suggests that the SED at A > 15 fiva is dom- 
inated by emission from this envelope. The SED model fits by 
Borges Fernandes et al. (2007) suggest that the system is seen 
under low to intermediate inclination (/ < 70°). 

The line profiles of the strong hydrogen recombination 
lines were modeled as emission from a spherically symmet- 
ric gas envelope (Benedettini et al. 1998). Acke et al. (2005) 
interpreted the profile of the optical [O i] line in the context of 
a wind originating from the surface layer of a passive disk. 

Using the AMBER visibilities measured at continuum 
wavelengths, we compare the ring diameter measured towards 
various PAs (covering a position angle range of ~ 68°) and 
find no indications for an elongation of the continuum-emitting 
region. This is consistent with the measured small Bry line 
closure phase (see Tab. 4), indicating that the brightness dis- 
tribution of the combined line- and continuum-emitting re- 
gion is nearly centro-symmetric, as in the case when the line- 
emitting region is seen nearly face-on. Furthermore, V921 Sco 
exhibits a continuum-emitting region which is more compact 
than expected for an irradiated dust disk, which might suggest 
that the NIR continuum emission is dominated by emission 
from an optically thick gaseous inner accretion disk (similar 
to MWC 297). 

Within the Bry line, we measure a slight increase in 
visibility (see Fig. 7). Furthermore, the continuum-corrected 
line-visibilities Vary (ranging between 0.7 and 0.5) show that 
the Bry-emitting region is also spatially resolved and only 
slightly more compact than the continuum-emitting region 
(7?Bry/^cont ~ 0.7). Therefore, the Bry region is too extended 
to be consistent with magnetospheric accretion or an X-wind 
as dominant Bry-emitting mechanism, which makes a strong 
stellar wind, a disk wind or a gaseous inner disk the most likely 
scenario. 
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Fig. 8. Top: The fitted ring radii for the continuum (Rcom. black points) 
and Bry hne (Rmy, red points) plotted as a function of stellar lumi- 
nosity. For the Bry line, we show the ring radius, determined by fit- 
ting only the spectral channel with the line center. For comparison, we 
plot the spatial extension of the stellar surface /?* (grey area) and the 
dust sublimation radius corresponding to dust sublimation tempera- 
tures of 2000, 1500, and 1000 K (computed using Eq. 2 and assuming 
£ = 1). Bottom: Providing a natural measure for the temperature dis- 
tribution within the circumstellar disk, we normalized Rg^ by the size 
of the continuum-emitting region Rcont- In particular, this normaliza- 
tion seems important considering the large range of stellar luminosi- 
ties covered by our sample. 



7. Discussion of general trends 

Since our sample covers a wide range of stellar parameters, 
we investigated whether the measured size of the continuum 
and Bry-emitting region correlates with the stellar parameters 
or the spectroscopic properties. Revealing such relationships 
could provide more insight into the involved physical mech- 
anisms and is also essential to confirm the empirically found 
correlations between the Bry luminosity and other estimators 
for the mass accretion rate (Calvet et al. 2004). In order to fur- 
ther expand our sample, we include in this section not only 
the five objects which we have investigated with AMBER, 
but also the MWC 480 Keck-Interferometer measurement pub- 
lished by Eisner (2007). For MWC 480, we assume L = 17 L©, 
d = 140 pc, Teff = 8700 K, a P-Cygni He-line profile (Acke 
et al. 2005), and a Bry line luminosity L(Bry)/Lo - -2.8 (as 
determined by L. Testi and A. Natta from unpublished TNG 
spectra; private communication). 
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7.1. Continuum emission 

Earlier interferometric investigations of the continuum- 
emitting region of HAeBe stars have established a relation be- 
tween the size of the continuum-emitting region and the stel- 
lar luminosity, suggesting that the continuum emission mainly 
traces hot dust located at the dust sublimation radius. In order 
to check for similar trends, we plot the determined A'-band con- 
tinuum ring radii as a function of the stellar luminosity L* in 
Fig. 8. For comparison, we computed the predicted inner disk 
radii corresponding to dust sublimation temperatures T^ubi of 
2000 K, 1500 K, and 1000 K using the analytic expression for 
the dust sublimation radius R^ubi by Monnier & Millan-Gabet 
(2002) 



^s, 



0.034 AU 



1500 K 



,subl 



Lq e' 



(2) 



In this relation, e denotes the ratio of the absorption efficien- 
cies of the dust at the sublimation temperature Tsubi and at the 
stellar effective temperature T*, i.e. e - Kp(T^ab\)lKp{Ti,) with 
Kp{T) = J Q.^^,^{A)B;^{T)dAI ^ B;i{T)dA. We assume e = 1, cor- 
responding to rather large dust grains of several /im in size and, 
therefore, to the innermost region of the inner dust rim. 

We find that the measured continuum radii of three of the 
five objects (HD 163296, HD 104237, HD 98922) follow the 

1 /2 

^subi ^ LJ. law rather closely assuming dust sublimation tem- 
peratures Teff between 1300 K and 1500 K, which is in agree- 
ment with the results obtained in earlier NIR broadband in- 
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terferometric survey observations (Monnier & Millan-Gabet 
2002). For the Herbig Be stars MWC297 and V921 Sco, the 
derived ring diameters are significantly more compact than the 
expected dust sublimation radii. A possible explanation for the 
small apparent continuum sizes might be inclination effects, 
although, based on complementary information from the lit- 
erature (see Sect. 6), we expect notable inclination only for 
HD 163296. Alternatively, the apparent disk size might be re- 
duced either due to gas absorption (allowing dust to exist closer 
to the star) or due to the emission of optically thick gas located 
inside of the dust-sublimation radius. As discussed in Sects. 6.4 
and 6.5, we favor these latter explanations for MWC297 and 
V921 Sco. 

7.2. Bry //ne emission 

7.2.1 . Correlations with stellar parameters 

Eisner (2007) speculated that the size of the Bry-emitting re- 
gion might depend on L*, with the (spatially compact) ac- 
cretion processes dominating for stars with low luminosity 
and (spatially extended) outflow processes dominating in the 
high-luminosity regime. Therefore, we investigate whether the 
derived physical size (RBry) or the relative size of the line- 
emitting region (RBry/Rcom, i-C- normalized by the size of the 
continuum-emitting region) scales with the stellar luminosity 
L*. As can be seen in Fig. 8, we do not find a clear correlation 
which includes all stars in our sample, although there might be 
a trend Rsry oc L^, if we exclude the Herbig Be star HD 98922. 
Likewise, there might be a relation between Rsry and the stel- 
lar temperature T* (Fig. 9, top), with a larger size of the Bry- 
emitting region for higher effective stellar temperature. The fact 
that we have to exclude HD 98922 seems to indicate that Bry 
traces not one unique physical mechanism for all stars in our 
sample, but at least two line-emitting mechanisms. As we argue 
in Sect. 6 and in the following subchapters, these mechanisms 
might include accretion and outflow processes, where the latter 
seem to be responsible for the tentative correlations with the 
stellar parameters noted above. However, given the small sta- 
tistical sample, more observations will be required to confirm 
these first trends. 

7.2.2. Correlations with spectroscopic parameters 

Concerning the spectroscopic parameters, we search for cor- 
relations between R^xy and RBry/Rmm and the Bry line lu- 
minosity L(Bry). Garcia Lopez et al. (2006) determined this 
quantity for all stars in our sample from the Bry line equiv- 
alent width of the circumstellar component. In principle, the 
accretion luminosity and mass accretion rate can be deter- 
mined from L(Bry) using the relation by Calvet et al. (2004) 
(logLacc/i-o = O.9x(log(LBry/i^0) + 4^-0.7)). However, since 
applying this relation for the B-type stars in our sample would 
require to extrapolate the empirical L(Biy)-L.^cc relation well 
outside the range of stellar properties (luminosity and mass) of 
the calibration sample, we refrain from this step and use the line 
luminosity instead. There seems to be a general trend, that with 
increasing line luminosity the size of the Bry-emitting region 



increases (Fig. 10). In 2001, Muzerolle et al. already proposed 
a correlation between W(Bry) (and thus also L(Bry)) and the 
size of the Bry-emitting region in the context of their mag- 
netospheric accretion model for T Tauri stars (the low-mass 
counterparts of HAeBe stars). In order to explain the empir- 
ical W(Bry)-Lacc correlation, they propose non-axisymmetric 
accretion, where the emitting area of discrete funnel flows in- 
creases with increasing accretion rate. However, since these 
funnel flows should still be located within a few stellar radii, 
we interpret the measurement of rather extended Bry-emitting 
regions in four out of six HAeBe stars as evidence against the 
hypothesis that Bry is predominantly tracing magnetospheric 
accretion in these stars. It is also interesting to note that the 
width of the Bry line seems to decrease with increasing stellar 
temperature (see Fig. 2). Stee & Bittar (2001) have shown that 
this could be explained by optical depth effects in wind models, 
although the influence of inclination efifects cannot be ruled out 
at the current time. 

Another spectroscopic parameter which was used exten- 
sively in the past to characterize the accretion and outflow pro- 
cesses around YSOs is the line profile of the Balmer lines, espe- 
cially Ha. Acke et al. (2005) classified the Ha line profile of all 
stars in our sample using the scheme "double-peaked", "single- 
peaked", "P-Cygni", and "inverse P-Cygni". The six stars in 
the considered sample exhibit three of these four Ha line pro- 
files (see Tab. 1), with the inverse P-Cygni profile being the 
exception^. As discussed in Sect. 1, interpreting these line pro- 
files without spatial information can lead to highly ambiguous 
results (for instance, Cidale & Ringuelet 1993 could already 
reproduce three of these four line profiles by varying the ve- 
locity law for an expanding wind). Thus, it seems promising to 
compare the spatially resolved information about the size of the 
Bry-emitting region with the shape of the Ha line profile, po- 
tentially revealing correlations between the spatial distribution 
and the kinematics of the emitting hydrogen. In this respect, it 
is important to notice that even under the assumption that Bry 
and Ha trace the same physical process, these lines will form 
in different spatial regions due to their different energy levels - 
with Bry tracing higher gas densities located closer to the star. 

As shown in Fig. 11, we find indications that RBry/Rmm 
depends on the Ha line profile classification: 
P-Cygni: For the HAeBe stars with a P-Cygni Ha-line profile 
(MWC 480, HD 98922), the Bry-emitting region is very com- 
pact (RBry/Rcont ^ 0.2, unrcsolvcd by the interferometer), con- 
sistent with an origin in a compact stellar wind, an X-wind or 
in magnetospheric accretion. P-Cygni Ha line profiles are gen- 
erally attributed to mass outflow (Finkenzeller & Mundt 1984; 
Catala & Kunasz 1987), although for T Tauri stars Muzerolle 
et al. (2001) showed that these profiles can also be reproduced 
in magnetospheric accretion models assuming blue-shifted ab- 
sorption by gas in an accretion-powered wind located outside 
of the magnetosphere. Therefore, the Bry line in MWC 480 
and HD 98922 might primarily trace magnetospheric accretion. 



^ As discussed, for example, by Walker (1972) and Sorelli et al. 
(1996), the inverse P-Cygni profile is generally believed to trace mass 
infall, but is found towards only ~ 5% of all HAeBe stars (Acke et al. 
2005) 
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while the optical Balmer lines are strongly affected by absorp- 
tion from outflowing gas. This interpretation is also consistent 
with the rather high accretion rates, which have been reported 
for these two objects (MWC480: logM^cc > -7.0 M0yr"', 
Eisner 2007; HD 98922: logM.^c = -5.76 Moyr^i, Garcia 
Lopez et al. 2006). 

Double-peaked: Several theoretical studies (e.g. Cidale & 
Ringuelet 1993) associate this Ha line profile with rotating 
gas flows, including disk winds. HD 163296, HD 104237, and 
V921 Sco, which fit into this classification, exhibit a Bry- 
emitting region which is only slightly more compact than 
the dust-sublimation radius (but still spatially resolved by 
our spectro-interferometric measurements). Considering the 
gaseous inner disk scenario, Muzerolle et al. (2004) modeled 
the continuum and line emission from such a disk and found 
that for low mass accretion rates (log Mace < -7 M0yr"'), the 
flux contribution from the gaseous inner disk is negligible. This 
leads us to favor the stellar wind or disk wind scenario for these 
three objects. 

Single-peaked: Exhibiting a single-peaked Ha-line profile, 
MWC297 shows a rather extended Bry-emitting region, ex- 
ceeding even the size of the continuum-emitting region. 
However, comparing the size of the Bry-emitting region around 
MWC297 to sources with double-peaked Ha-line profiles 
(such as V921 Sco), we do not find fundamental differences 
within the measurement uncertainties (especially taking the 
indications for a particularly compact continuum-emitting re- 
gion around MWC 297 into account, as discussed in Sect. 6.4). 
Theoretical models have associated single-peaked Ha line pro- 
files with stellar winds or with disk winds (Cidale & Ringuelet 
1993), especially for massive YSOs (Sim et al. 2005). 

Based on the currently available spectroscopic and long- 
baseline spectro-interferometric data, it seems not feasible to 
discern between stellar wind and disk wind scenarios. This is 
due to the fact that the gas velocity field in stellar winds as 
well as in disk winds, comprise both rotating and expanding 
velocity components, which can result in similar line profiles 
as well as in similar spatial scales for the line-emitting region. 
To solve these ambiguities, future observations will require a 
very good sampling of the line visibility function in order to 
measure the subtle imprints resulting from the differences of 
the gas kinematics close to the star. 

7.2.3. Implications on the empirical L(Br7)-Lacc 

relation and the accretion-outflow connection 

Within the last decade, a correlation between the Bry line lu- 
minosity and the mass accretion rate (as independently de- 
termined from the UV excess, which is commonly attributed 
to the accretion shock) could be established for pre-main- 
sequence stars with masses between ~ 0.01 (Natta et al. 2004) 
and ~ 4 Mq (Muzerolle et al. 1998b; van den Ancker 2005). 
Being applicable even in regions of high extinction and for 
very low mass accretion rates, this correlation is of consid- 
erable practical importance, although neither the traced pro- 
cesses nor the underlying physical mechanisms are yet iden- 
tified. It was suggested that the Bry line emission might be 



a direct tracer of mass accretion, and that for large accretion 
rates (M-^cc > 10"^ MQyr"') the contributions from other pro- 
cesses might be negligible (van den Ancker 2005). However, 
attempts to reproduce the L(Bry)-Lacc relation using magneto- 
spheric accretion models were not successful, but found instead 
that the Bry line luminosity depends mainly on other parame- 
ters such as the size of the magnetosphere and the gas temper- 
atures (Muzerolle et al. 1998a,b). 

Based on our spatially resolved observations of the Bry- 
emitting region around five Herbig Ae/Be stars, we cannot sup- 
port the scenario that Bry generally traces gas located in the 
magnetospheric infall zone. Considering only those stars in our 
sample which are within the current calibration range of the re- 
lation (i.e. rejecting the early Herbig Be stars), we find that 
for one star (HD 98922) the size of the Bry-emitting region 
is consistent with magnetospheric accretion, while two stars 
(HD 104237, HD 163296) exhibit a Bry-emitting region which 
is more consistent with an extended stellar wind or a disk wind 
scenario (/^Bry ^ 0.6 Rcom)- Therefore, our observations im- 
ply that, at least for some HAeBe stars, Bry is not a primary 
tracer of accretion, but likely indirectly linked to the accretion 
rate, e.g. via accretion-driven mass loss. This ambivalent ori- 
gin of the line emission (tracing partially infalling and outflow- 
ing matter) makes the empirical correlation between L(Bry) 
and Lace even more remarkable and suggests a tight, quanti- 
tative connection between the accretion and ejection processes 
in YSOs. The existence of such a connection has already been 
suggested based on measurements of the luminosity of accre- 
tion and outflow-tracer lines, yielding an ejection efficiency of 
Mejeet/Maee ~ 0. 1 -0.2 (Cabrit ct al. 1990; Hartigan ct al. 1994). 

7.2.4. Implications on the jet-launching mechanism for 
Herbig Ae/Be stars 

Our measurement of rather extended Bry-emitting regions also 
has important implications on the launching mechanism for 
HAeBe jets. Although most stars in our sample show only indi- 
rect signatures of outflow activity (e.g. CO outflows, outflow- 
tracing line emission), HD 104237 and HD 163296 are also 
known to be associated with coUimated microjets. On large 
scales, the interaction region between the jet ejecta and ambient 
material is traced by Herbig-Haro objects, while on arcsecond- 
scales, the jets were also imaged with the HST in forbidden 
lines and in the hydrogen Ha and Lya transitions (Grady et al. 
2004; Wassell et al. 2006). 

Therefore, it seems likely that our AMBER Bry obser- 
vations trace ionized gas located at the base of these micro- 
jets, just experiencing the acceleration and magnetic confine- 
ment processes which are required to achieve the highly colli- 
mated gas flows detected by HST. Since it is expected that stel- 
lar winds, lacking any collimation mechanism, result in wide- 
angle outflows, we clearly favor the disk wind model for these 
objects in order to explain the observational evidence for (a) a 
rather extended Bry-emitting region on sub-AU scales and, si- 
multaneously, (b) collimated narrow-angle microjets on scales 
of several hundred AU. This main conclusion is in qualitative 
agreement with HST spectroscopic observations, which mea- 
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sured the jet opening angle and the toroidal velocity field on 
spatial scales ~ 50 AU from the base. Traced back to the foot- 
point of the jet, the measured velocity fields seem more consis- 
tent with a launching of the jet from the disk and not an X-wind 
or stellar wind (Bacciotti et al. 2002; Testi et al. 2002; CoflFey 
et al. 2007). 

Future spectro-interferometric observations will allow us to 
confirm this hypothesis by measuring the spatial distribution 
and kinematics of the line-emitting gas and comparing it with 
the orientation and velocity of the large-scale outflow structure. 



8. Conclusions 

We summarize the results from our investigation using 
VLTl/ AMBER data of five HAeBe stars as follows: 

1. The /T-band continuum-emitting region was spatially re- 
solved for all objects. For the Herbig Ae and late Herbig Be 
stars the characteristic size of this region scales roughly 
with the square-root of the stellar luminosity, supporting 
the scenario that for these sources, the continuum is domi- 
nated by thermal emission from the inner rim of the dust 
disk. For the early-type Herbig Be stars MWC297 and 
V921 Sco, on the other hand, the continuum-emitting re- 
gion is significantly more compact than predicted by the 

1 11 

RavM ^ l^i, relation, maybe suggesting that for these ob- 
jects the NIR continuum is dominated by hot gas emission. 
Fits of inclined ring geometries to the visibilities measured 
on HD 163296 indicate object elongation, consistent with 
a notably inclined disk geometry, while for HD 98922 and 
V921 Sco we find no indications for a large disk inclination. 

2. Closure phases were measured for four objects in the K- 
band continuum regime. For V921 Sco we measured the CP 
also within the spectrally resolved Bry line, suggesting a 
nearly centro-symmetric line-emitting region (in agreement 
with the low inclination angle). 

3. In the spectrally resolved Bry lines of HD 163296, 
HD 104237, HD 98922, and V921 Sco, we measure an in- 
crease of visibility, indicating that the line-emitting re- 
gion is more compact than the continuum-emitting region. 
MWC 297 shows a decrease of visibility within the Bry- 
line, although this is likely related to the compactness of 
the continuum-emitting region and not to fundamental dif- 
ferences in the line formation mechanism. 

4. Judging on the currently available limited statistics, we do 
not find a trend which relates the size of the Bry line- 
emitting region with the stellar luminosity, as was pro- 
posed by Eisner (2007). Instead, the continuum-normalized 
size of the line-emitting region Rsry/Rcont might correlate 
with the stellar effective temperature, although more obser- 
vations will be required to confirm this trend. Maybe the 
most remarkable correlation we found was that stars with 
a P-Cygni Ha line profile and a high mass-accretion rate 
seem to show particularly compact Bry-emitting regions 
(^Bry/^cont < 0.2), whilc stars with a double-peaked or 
single-peaked He-line profile show a significantly more ex- 
tended Bry-emitting region (0.6 < Rsry /Rcont i 1-4). 



5. Only for HD 98922, the emitting-region is compact enough 
to support the hypothesis that most of the Bry-emission 
emerges from magnetospheric accretion columns. In the 
other cases, the Bry-emitting region is only slightly more 
compact than the dust sublimation radius, supporting the 
idea that the emission mainly emerges from an extended 
stellar wind or a disk-wind. This makes it unlikely that Bry 
primarily traces magnetospheric accretion in HAeBe stars, 
suggesting that even for Herbig Ae stars, this line is only an 
indirect tracer of the mass accretion rate. 

Acknowledgements. We are grateful to A. Chelli and G. Duvert for 
carrying out independent checks on the AMBER data reduction. 
Furthermore, we would like to thank the referee for his comments, 
which helped to improve the paper. This publication makes use of 
data products from the Two Micron All Sky Survey, which is a joint 
project of the University of Massachusetts and the Infrared Processing 
and Analysis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and the National 
Science Foundation. NSO/Kitt Peak FTS data used here were pro- 
duced by NSF/NOAO. 

References 

Acke, B., van den Ancker, M. E., & Dullemond, C. R 2005, 

A&A, 436, 209 
Acke, B., Verhoelst, T, van den Ancker, M. E., et al. 2008, 

A&A, 485, 209 
Bacciotti, F, Ray, T. P, Mundt, R., Eisloffel, J., & Solf, J. 2002, 

ApJ, 576, 222 
Benedettini, M., Nisini, B., Giannini, T., et al. 1998, A&A, 339, 

159 
Blandford, R. D. & Payne, D. G. 1982, MNRAS, 199, 883 
Blondel, P F C. & Djie, H. R. E. T. A. 2006, A&A, 456, 1045 
Bohm, T., Catala, C, Balona, L., & Carter, B. 2004, A&A, 427, 

907 
Borges Fernandes, M., Kraus, M., Lorenz Martins, S., & de 

Araujo, E X. 2007, MNRAS, 377, 1343 
Brittain, S. D., Simon, T., Najita, J. R., & Rettig, T W. 2007, 

ApJ, 659, 685 
Cabrit, S., Edwards, S., Strom, S. E., & Strom, K. M. 1990, 

ApJ, 354, 687 
Calvet, N., Muzerolle, J., Briceno, C, et al. 2004, AJ, 128, 

1294 
Catala, C, Donati, J. F, Bohm, T., et al. 1999, A&A, 345, 884 
Catala, C. & Kunasz, P B. 1987, A&A, 174, 158 
Causi, G. L., Antoniucci, S., & Tatulli, E. 2007, ArXiv e-prints, 

711 
Cidale, L. S. & Ringuelet, A. E. 1993, ApJ, 411, 874 
Coff^ey, D., Bacciotti, F, Ray, T P, Eislofi'el, J., & Woitas, J. 

2007, ApJ, 663, 350 
Corcoran, M. & Ray, T P 1998, A&A, 331, 147 
de Winter, D. & The, P S. 1990, Ap&SS, 166, 99 
Donati, J.-F, Semel, M., Carter, B. D., Rees, D. E., & Collier 

Cameron, A. 1997, MNRAS, 291, 658 
Drew, J. E., Busfield, G., Hoare, M. G., et al. 1997, MNRAS, 

286,538 
Dullemond, C. R, Dominik, C, & Natta, A. 2001, ApJ, 560, 

957 



Kraus et al.: Origin of hydrogen line emission in Herbig Ae/Be stars 



15 



Eisner, J. A. 2007, Nature, 447, 562 

Eisner, J. A., Lane, B. R, Hillenbrand, L. A., Akeson, R. L., & 

Sargent, A. I. 2004, ApJ, 613, 1049 
Ferreira, J. 1997, A&A, 319, 340 

Ferreira, J., Dougados, C, & Cabrit, S. 2006, A&A, 453, 785 
Finkenzeller, U. 1985, A&A, 151, 340 
Finkenzeller, U. & Mundt, R. 1984, A&AS, 55, 109 
Garcia Lopez, R., Natta, A., Testi, L., & Habart, E. 2006, A&A, 

459, 837 
Grady, C. A., Devine, D., Woodgate, B., et al. 2000, ApJ, 544, 

895 
Grady, C. A., Perez, M. R., Talavera, A., et al. 1996, A&AS, 

120, 157 
Grady, C. A., Woodgate, B., Torres, C. A. O., et al. 2004, ApJ, 

608, 809 
Habart, E., Testi, L., Natta, A., & Vanzi, L. 2003, A&A, 400, 

575 
Hartigan, R, Morse, J. A., & Raymond, J. 1994, ApJ, 436, 125 
Henning, T., Burkert, A., Launhardt, R., Leinert, C, & 

Stecklum, B. 1998, A&A, 336, 565 
Higdon, S. J. U., Devost, D., Higdon, J. L., et al. 2004, R\SP, 

116,975 
Hutsemekers, D. & van Drom, E. 1990, A&A, 238, 134 
Isella, A. & Natta, A. 2005, A&A, 438, 899 
Isella, A., Testi, L., Natta, A., et al. 2007, A&A, 469, 213 
Kraus, S., Preibisch, T., & Ohnaka, K. 2008, ApJ, 676, 490 
Leinert, C, van Boekel, R., Waters, L. B. R M., et al. 2004, 

A&A, 423, 537 
Malbet, R, Benisty, M., de Wit, W.-J., et al. 2007, A&A, 464, 

43 
Malfait, K., Bogaert, E., & Waelkens, C. 1998, A&A, 331, 211 
Mannings, V. & Sargent, A. L 1997, ApJ, 490, 792 
Mathis, J. S. 1990, ARA&A, 28, 37 
Mestel, L. 1968, MNRAS, 138, 359 
Millan-Gabet, R., Schloerb, R P, & Traub, W. A. 2001, ApJ, 

546, 358 
Monnier, J. D., Berger, J.-P, Millan-Gabet, R., et al. 2006, ApJ, 

647, 444 
Monnier, J. D. & Millan-Gabet, R. 2002, ApJ, 579, 694 
Monnier, J. D., Millan-Gabet, R., Billmeier, R., et al. 2005, 

ApJ, 624, 832 
Muzerolle, J., Calvet, N., & Hartmann, L. 1998a, ApJ, 492, 743 
MuzeroUe, J., Calvet, N., & Hartmann, L. 2001, ApJ, 550, 944 
Muzerolle, J., D'Alessio, P., Calvet, N., & Hartmann, L. 2004, 

ApJ, 617,406 
Muzerolle, J., Hartmann, L., & Calvet, N. 1998b, AJ, 1 16, 2965 
Natta, A., Palla, R, Butner, H. M., Evans, II, N. J., & Harvey, 

R M. 1993, ApJ, 406, 674 
Natta, A., Prusti, T., Neri, R., et al. 2001, A&A, 371, 186 
Natta, A., Testi, L., Muzerolle, J., et al. 2004, A&A, 424, 603 
Nisini, B., Milillo, A., Saraceno, P, & Vitali, R 1995, A&A, 

302, 169 
Oudmaijer, R. D. & Drew, J. E. 1999, MNRAS, 305, 166 
Petrov, R. G., Malbet, R, Weigelt, G., et al. 2007, A&A, 464, 1 
Pogodin, M. A. 1994, A&A, 282, 141 
Pudritz, R. E. & Norman, C. A. 1983, ApJ, 274, 677 
Richichi, A., Percheron, I., & Khristoforova, M. 2005, A&A, 

431,773 



Shu, R, Najita, J., Ostriker, E., et al. 1994, ApJ, 429, 781 
Sim, S. A., Drew, J. E., & Long, K. S. 2005, MNRAS, 363, 615 
Sorelli, C, Grinin, V. P, & Natta, A. 1996, A&A, 309, 155 
Stee, P & Bittar, J. 2001, A&A, 367, 532 
Strafella, R, Pezzuto, S., Corciulo, G. G., Bianchini, A., & 

Vittone, A. A. 1998, ApJ, 505, 299 
Tannirkulam, A., Monnier, J. D., Millan-Gabet, R., et al. 2008, 

ApJ,677,L51 
Tatulli, E., Isella, A., Natta, A., et al. 2007a, A&A, 464, 55 
TatuIIi, E., Millour, R, Chelli, A., et al. 2007b, A&A, 464, 29 
Testi, L., Bacciotti, R, Sargent, A. I., Ray, T. P, & Eisloffel, J. 

2002, A&A, 394, L31 
van Boekel, R., Waters, L. B. R M., Dominik, C, et al. 2003, 

A&A, 400, L21 
van den Ancker, M. E. 2005, in High Resolution 

Infrared Spectroscopy in Astronomy, ed. H. U. Kaufl, 

R. Siebenmorgen, & A. R M. Moorwood, 309—1- 
van den Ancker, M. E., Bouwman, J., Wesselius, P. R., et al. 

2000, A&A, 357, 325 
van den Ancker, M. E., de Winter, D., & Tjin A Djie, H. R. E. 

1998, A&A, 330, 145 
van den Bergh, S. & Herbst, W. 1975, AJ, 80, 208 
Vinkovic, D. & Jurkic, T. 2007, ApJ, 658, 462 
Walker, M. R 1972, ApJ, 175, 89 
Wang, S. & Looney, L. W. 2007, ApJ, 659, 1360 
Wassell, E. J., Grady, C. A., Woodgate, B., Kimble, R. A., & 

Bruhweiler, R C. 2006, ApJ, 650, 985 
Weigelt, G., Kraus, S., Driebe, T., et al. 2007, A&A, 464, 87 



List of Objects 



'HD 101328' on page 3 
'HD 118934' on page 3 
'HD 135382' on page 3 
'HD 159941' on page 3 
'HD 171960' on page 3 
'HD 177756' on page 3 
'HD 104237' on page 4 
'HD 163296' on page 4 
'HD 98922' on page 4 
'MWC 297' on page 4 
'V921 Sco' on page 4 



16 



Kraus et al.: Origin of hydrogen line emission in Herbig Ae/Be stars 



a) V{A) 
velocity [km/s] 

-1000 1000 




> 

E- 



1 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 

? 

0.9 

0.8 

0.7 

0.6 

i 0.5 

' 0.4 

^ 0.3 

5 0.2 

0.1 


9 

¥ ^ 

E 7 






CO 

Q 



PEL 35.1 m,PA7r, R=500 

■ PBL58.0m,PA78°, R=1500 

PEL 58.5 m, PA 84°, R=500 



PBL35.1 m,PA71°, R-500 
PEL 58.0 m, PA 78°, R-1500 
PEL 58.5 m, PA 84°, R=500 




PBL35.1 m,PA7l'', R=500 
PEL 58.0 m, PA 78°, R=1500 
PBL 58.5 m, PA 84°, R-500 



— 1_ 



— 1_ 



PBL 35.1 m, PA 71 °, R=500, continuum 
-B-- PBL 58.0 m, PA 78°, R=1500, continuum 

- "- PBL 58.5 m, PA 84°, R=500, continuum 
— ^ PBL 35.1 m, PA 71 °, R-500, line ; 
-B-- PBL 58.0 m, PA 78°, R-1500, line ; 

- ■* - PEL 58.5 m, PA 84°, R-500, line 



;?lgJ5J5^p^?n;^j|j^|{] jt 



i)^' 



2.16 2.17 

Wavelength [urn] 





1 

■ 0.8 
0.6 
0.4 
0.2 




1e-09 


b) SED 


Kumcz Atmosphere A4 star (8000K, 9=4.04) 

DereddenedSED(A^.3l'^;Rv=3.r) -^ 


1e-10 


/v. ^ 


1e-11 


^ -^^r---. ■ 


1e-12 

l6-t3 


\ 

\ 

\ 
\ 
\ 
\ 



1 10 

A[Mm] 

c) Vcont(B), continuum 

Spadai Frequency [10^ K] 
10 20 30 40 50 

"".,„, ' ' IJD 9=7.47 mas 

'■■';", Gaussian 6=4.93 mas 

''■-::■. Ring D=4.97 mas 

'v.. 1 AMBER data, continuum around Bry 



100 120 



Prcijected Baseline Length [m] 

d) VDr7(B), line 

Spadai Frequency [10^ k] 















1 Ring D=3:i4 mas 
■ .. f Bryline -•- J 


\ 













100 120 



Projected Baseline Length [r 



e) Geometric Models 

dRA [AU] 
0.4 0.2 0.2 0.4 




54321012345 
dRA [mas] 



Fig. 3. Spectroscopic, spectro-interferometric, and photometric data for HD 104237: In (a), we show the Bry spectra extracted from our ISAAC 
and AMBER data (top panel), the measured wavelength-dependent visibilities Vtoi for different baselines (2nd panel, the data points are shown 
with statistical errors bars, whereas the estimated calibration errors are shown in the bottom-right corner), the continuum-corrected Bry-line 
visibility V^^^ (3rd panel), the UD diameter ©tm derived from Vto, (4th panel), and the diameters £)cont and D^^y derived for the continuum- and 
the Bry-emitting region using a 2-ring model (5th panel). In (b), the SED is shown, including photometric data from the literature (grey points), 
archival ISO spectra (magneta), and archival 5p;fzer/IRS spectra (blue). In (c), we plot the measured A'-band continuum-visibility as a function 
of baseline length and in (d) the continuum-corrected Bry-line visibility. Panel (e) shows the derived ring diameters for the continuum- and 
line-emitting region plotted as function of position angle. 
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Fig. 4. Spectroscopic, spectro-interferometric, and photometric data for HD 163296 (similar to Fig. 3). 
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Fig. 5. Spectroscopic, spectro-interferometric, and photometric data for HD 98922 (similar to Fig. 3). 



> 



m 
> 



Kraus et al.: Origin of hydrogen line emission in Herbig Ae/Be stars 

a) V{A) . .. b) SED 



19 



velocity [km/s] 

-1000 



3.5 
V 3 
£ 2.5 
lT 2 

1.5 

1 

1 

0.9 

0.8 

S 0.7 

S 0.6 

5 0.5 

g 0.4 

-* 0.3 

0.2 

0.1 





5 ■ 
4 ■ 
3 
2 
1 



Q 



7 ■ 
6 ■ 
5 
4 

3 ■ 
2 ■ 
1 



AMBER R-1500 

■ ISAAC R-1500 

■ ISAAC R^8900 



H, PA 40 , R-1500 



PEL 43.8 m, PA 40°, R=1500 



n, PA 40°, R-1500 



-U-^iii-^»,i-*ft + ,-»i-i+ 



PBL43.8 m, PA 40°, R-1500, continuum 
PEL 43.8 m, PA 40°, R-1500, line 



-»U-tfi+l-M-i+l4i-l* 



3 




# 


0. 


2.5 








? 








1.5 
1 


=3 
< 

i 






0.5 






5 

4 









3 


2 






2 


1.5 


5" 
< 

s 

Q 


E, 

O 

LU 
□ 


1 


1 



Kurucz Atmosphere (T=24000K, g=4.04) 
Deteddened SED (Av=8'^; Rv=3.l") 




c) Vcont(B), continuum 

Spadal Frequency [10^ K] 



UD 9=7.96 mas ■ 
Gaussian ft=5.07 mas ■ 
Ring D=5.40 mas 
AMBER data, continuum around Bry - 



Prcijected Baseline Length [m] 

d) VDr7(B), line 

Spatial Frequency [10^ k] 









Ring D=7.55 mas 

Bry line — •- 






\ 








\ 



Projected Baseline Lengtti [m] 

e) Geometric Models 

dRA [AU] 
0.5 0.5 1 




54321012345 



Wavelength [urn] dRA [mas] 

Fig. 6. Spectroscopic, spectro-interferometric, and photometric data for MWC 297 (similar to Fig. 3). 
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Fig. 7. Spectroscopic, spectro-interferometric, and photometric data for V921 Sco (similar to Fig. 3). 



